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a^ ■ ABSTRACT 



An exactly solvable model of a quantum spin interacting with a spin environment 
is considered. The interaction is chosen to be such that the state of the environment is 



O^' conserved. The reduced density matrix of the spin is calculated for arbitrary coupling 



strength and arbitrary time. The stationary state of the spin is obtained explicitely in 



c^ ■ the t ^ oo limit. 



The problem of a quantum system interacting with a heat bath has been extensively 
studied in various contexts during the last few decades. Originating from quantum 
optics in connection with studying spontanious emission and resonant adsorpsion [1-2], 
it has become an important issue in condensed matter physics. The first, and probably 
one of the most famous works in this field was by Caldeira and Leggett who studied 
effects of dissipation on the probability of quantum tunneling [3] . In this problem the 
heat bath is modeled by a set of noninteracting harmonic oscillators linearly coupled 
to the quantum system. This model of heat bath has been mathematically justified 
in [4], and has become a widely accepted description of dissipative quantum dynamics, 
which advantageously combines both microscopical and phenomenological aspects of 
interaction between a quantum system and phonons or delocalized electrons [5]. A 
similar problem of a system in a double well potential under the influence of a Heat 
Bath was studied in connection with magnetic flux tunneling in Josephson junctions [6]. 
It has been shown that the system looses its quantum coherence due to the interaction 
with the heat bath and, in case of zero temperature and sufficiently strong coupling, 
completely localizes in one of the wells. Later this model was formalized by the spin- 
boson Hamiltonian [7,8] that received a lot of attention in modern condenced matter 
literature. 

The study of quantum dissipation effects in the case of coupling to the fermionic 
heat bath has also received some attention in the literature [9]. For Hubbard-like cou- 
pling, behaviour quite different from that of the bosonic heat bath, has been found. The 
case of magnetic coupling similar to RKKY interactions [10] was also extensively studied 
in conjecture with magnetic grains and giant spins of macromolecules interacting with a 
spin environment [11-13]. It is believed that this type of heat bath can not be mapped 
onto a bosonic bath model and needs separate treatment [11,12]. This mechanism of 
quantum relaxation turns out to be effective especially at low temperatures resulting 



in a set of interesting phenomena, such as "degeneracy blocking" caused by the nuclear 
spins [11]. 

Another interesting and quite general effect that results in the ineraction of quan- 
tum system with its environment is the destruction of quantum interference in the 
quantum system due to such interaction. This process, usually termed in the literature 
as decoherence, has attracted attention of both theorists and experimentalists not only 
due to its fundamental importance in quantum mechanics, but also due to the new fastly 
developing fields, such as quantum computing and quantum information theory, where 
decoherence is one of the major obstacles on the way of practical realizations of various, 
presently mostly hypothetical devices, such as quantum computers, etc [14]. No matter 
how well such a device is isolated from its environment, being essentially a macroscopic 
system, it will inevitably interact with the environment, resulting in the loss of interfer- 
ence between the states and thus desrupting its proper functioning. Several models have 
been proposed to study properties of decoherence in quantum computers and in more 
general systems [14-16]. The essential feature of these models is that the interaction is 
set up in such a way that there is no energy exchange between the system under con- 
sideration and the environment, so the system's energy is conserved. This corresponds 
to a situation when a quantum system is very well isolated from its evironment, so only 
a phase exchange is allowed. One considers a hamiltonian 



H = Hs + Hb+V, (1) 

where the first term Hs in (1) corresponds to the system alone, Hb is the heat bath and 
V is the interaction between the bath and the system. It is assumed that the expectation 
value of Hs is conserved during the system's evolution. This can be formalized by the 
assumption that the system's hamiltonian Hs commutes with the full hamiltonian in 



(1), and in particular with the interaction term V. Hs and Hb naturally commute with 
each other as they act in different subspaces. This feature often allows one to carry out 
exact solutions for the system's reduced density matrix (to be defined later) in the basis 
of eigenavalues of the self-hamiltonain of the system Hs ■ 

In this work we consider an opposite extreeme, a model having a property that 
the state of the bath is preserved. The motivation to study such a model comes from 
a very common assumption in the literature on the propeties of the heat bath. It is 
often assumed that the heat bath has so many degrees of freedom that the effects of the 
interaction with the system dissipate away in it and will not infiuence the system back 
to any significant extent so that the bath remains described by a thermal equilibrium 
distribution at constant temperature, irrespective of the amount of energy and polar- 
ization diffusing into it from the system [1,2]. This assumption is also called mollecular 
chaos or Stosszahlansatz. There are some attempts to analyze this assumtion [17], but 
we beleive it is still little understood. This work is aiming to contribute to this topic. 

Usually the following picture is assumed: the initial state of the full system is given 

by 



p(0)=Ps(0)®pb(0), (2) 

where ps(0) and Pb(0) are the initial density matricies of the system and the bath 
respectively. The heat bath is initially assumed to be in thermal equilibrium state and 
the two systems are not initially entangled. When the interaction is switched on at time 
t = 0, the full system's evolution is given approximately by the following 



p{t) - ps{t) ® pb{0) , (3) 
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that is, the state of the heat bath does not change in time to any significant extent. 
This assertion, formahzed by the Markoffian approximation, leads to the famous Pauh 
master equations, which provide a key to understanding many profound phenomena in 
quantum optics and condensed matter physics [1,2]. 

Here we propose a model for which equation [3] is exact. Note that the form of 
interaction V between the system and the bath has not been specified yet. Comparing 
equations (1) and (3), one can notice that in order for (3) to hold, one can simply require 
the following commutation relation 

[V,Hb] = [H,Hb]=0. (4) 

Let us comment on this assumption. The most common form of the interaction 
between the system and the bath is the coupling between operators Qi, acting in the 
subspace of the system and the bath operators F^, so that 

V = J2Qr®F,. (5) 

i 

The choice of operators Qi and Fi is usually determined by the particular features 
of the physical situation under consideration. However usually operators Qi and Fi 
are not diagonal in the energy representation of the system and the bath and so the 
relation (4) may not be rigorously satisfied in general and is usually postulated by the 
Stosszahlansatz assumption. 

In this paper we study a model for which relation (4) is satisfied directly due to 
the commutation properties of the interaction and bath hamiltonians. Such a property, 
even though not very common in the literature, seems worth exploring due to the above 
arguments. Moreover, this model allows exact solution for the reduced density matrix. 



which is a rather rare example in the hterature on decoherence and thermahzation of 
quantum systems. 

In particular, we consider a model of interaction between a two-level quantum 
system and spin environment [11-14]. A bath of noninteracting spins a^ is coupled to 
the two-level system a^ under consideration. The coupling is chosen to be such that the 
Hamiltonian for the full system (two- level system -|- spin bath) is: 

H = AaO + Y. u^a\ + a^ J^ g^a^ , (6) 

k k 

The first term in (6) corresponds to the two-level system and we will refer to it in 
the following as the central spin; the second is the spin bath, and the last one is the 
interaction between the bath and the two-level system. Here 2A is the bare magnetic 
resonance (MR) frequency of the central spin, Uk and gk are the frequencies and the 
coupling constants respectively for the spins of the spin bath. The spin bath self- 
hamiltonian obviosly commutes with the full hamiltonian in (6) and so the state of the 
bath is conserved in the course of the full system's evolution. 

The interaction is assumed to be switched on at time t = and the two systems 
(the central spin and the bath) are initially not entangled with each other. The density 
matrix of the full system is given by 

p(0)==(|lo)(lo|)®^e-/3^-. (7) 

Here P = l/ksT is the inverse temperature and by Hb we denote the self-hamiltonian 
of the bath, which is the second term in (6) . The spin bath is assumed to be initially 
in thermal equilibrium at temperature T. Z is the normalization constant or partition 
function of the free heat bath, i.e., a system of noninteracting spins 
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Z = Tt e"^^^ = n (2 cosh (3uJk) ■ (8) 

k 

To avoid unnecesary mathematical complications we have assumed that the spin is 
initially in the excited state, even though the calculation can be carried out exactly for 
arbitrary superposition of both ground and excited states. The full system evolves in 
time t quantum mechanically according to 

p{t)=U{t)p{0)U-\t). (9) 

Here and in the following we assume that h = 1. 

The evolution operator U{t) = e~'^^^ can be explicitely calculated by expanding 
the exponent and combining the terms, that correspond to the expansion of cosine and 
sine respectively, resulting in 



U = e- 



-HbI 



^ .AaO + OaS . ^ 
cos 'ft — t sm 71 
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(10) 



where 



Q^J29k^'^ (11) 



7=(A2 + 02)2_ (12) 

The reduced density matrix for the central spin is given by 

p^{t)=Tr'[p{t)], (13) 

where prime denotes that the trace is taken over the states of the spin bath only. 



We calculate the elements of the density matrix using the following technique. 
Consider the magnetization of the central spin, which is given by the expectation value 
of a^ or the difference of the reduced density matrix diagonal elements 



(^°W) = (lo|p'^(t)|lo)-(Oo|p'^(t)|Oo) 



(14) 



With the use of (7), (9)-(10), equation (14), after some algebra, can be rewritten as 



(^°W) = ^ 



-PHe 



cos^ 7t 



02 -A^ 



sin 7t 



(15) 



z V r 

The trace in the basis of eigenstates of Hb and O becomes a sum over Ising-like variables 

Sk = ±1 



(-;w) = E^A,^"^'' 



{sfe} 



V 1 
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(16) 



where the sum is taken over all possible configuartions of Sk- In equation (14) by Hb 
and O (see eqs.(ll)-(12)) we mean, of course, the eigenvalues of these operators, i.e., 
Hb = Ylik ^kSk and O = ^^ Qk^k- In (16) we have also intoduced operator A^ given by 



A. 



d_ 

drj 



/•2t 

+ AM di], 

r]=2t Jo 



(17) 



which can be conviniently interchanged with the summation in (16). In order to compute 
the sum (16) we employ the following identity [18]: 



Re f e^"" Jo (A^/7]^-xA dx 



sin?7(A2 + 02y 



(18) 



(A2 + 02)2 

Here Jo{z) is the zeroth order Bessel function. At this point performing the summation 
in (16) is a straightforward procedure as O enters linearly in the exponent in (18) and 



the sum is equivalent to the simple calculation of the partition function for a system of 
uncoupled spins in external magnetic field. So using (18) and also (11)-(12), after some 
algebra we obtain: 



(a°(t)) = A, Re T dxJ^ {^^/^f^^) | Y. 

{Sfej 

Ar, Re / dxJo (a^t]^ - xA ^x 



-pHs+ixVL 



and 



$(x) = TT [cos (7jt.a; — itanh/Jwfc sin 



9kX\ 



(19) 



(20) 



The trick used in (18)- (20) is simular to that in calculation of partition functions for 
certain mean field models, such as Curie- Weiss model; see [19] and references therein. 

The off-diagonal elements of the reduced density matrix can be calculated in a 
simular way. The final result is 



p\,{t)=K'lra 



$(a;) 



A I X 

'o 



V Ji i^i 



1]^ — X' 



rj^ — x^ 



$(a;) dx 



(21) 



where 



A' 



p2t 

+ A dr]. 

ri=2t Jo 



(22) 



Equations (19)-(21) constitute the main result of this work. All calculations up to this 
point were exact for arbitrary coupling constants and energy splittings of the spins. 

Let us specify the coupling constants gk and MR frequencies 2uJk of the external 
spins a^. For simplicity we assume that the energy splittings for all "external" spins 
are the same Uk = uj and the coupling constants are randomly distributed with average 
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(g) and second moment {g'^)- Expression (20) for $(x) can be exponentiated thus 
transforming the product into summation in the exponent 



$(a;) = exp 



5Z^fc(^) 



(23) 



where 



Ak{x) = 7T In [cos^ gkX + tanh /3uj sin^ gkx] — itan ^ [tanh /3uj tangkx] . (24) 

Assume now that the sum in (23) contains A^ terms, that is, the spin bath consists of 
N spins, where N ^ 1. Ak{x) are random numbers with some average {A{x)). with 
this assumption expression (23) becomes 



^{x) =expA(A(x)) 



(25) 



where ( ) denote the average taken over the couphng constants. In this work we assume 
that {g) = and {g"^) = ^, where C is of order unity. This choice is made only for 
simphcity of calculations and obviosly other possibilites for distributions of uJk and gk 
can be explored. 

Averaging of (24) can be done easily by expanding it up to the second order in g 
with the result 



$(a;) = exp 



C 



-X 



(26) 



2 cosh (3uj 

This relation, when inserted into (19)-(21), gives the explicit analytical expressions for 
the density matrix of the spin. The magnetization (the difference between the diagonal 
elements) is an oscillatory function, which decays to the limiting value for (^"^(oo)). 
This can be explicitly calculated and after the straightforward manipulations such as 
change of order of integration in (19) one obtains: 
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((tO(oo)) = Acosh/3a;-^/^exp (z^) erfc{z) (27) 



where 



z = ^=cosh/3u;. (29) 

V2C 

Here erfc{z) is the comphmentary error function. A similar straightforward calculation 
shows that the off-diagonal elements vanish for t = oo. 

One could expect that the central spin "thermalizes" due to the interaction with 
the spin bath, i.e., its density matrix reaches the state distributed according to the 
Boltzmann distribution. One can see that for this model this is not the case. Moreover, 
it is a straightforward observation that p^(oo) depends on the initial state p^(0) and 
thus the conservation of the state of the bath is not sufficient to represent molecular 
chaos assumption and cannot represent a "true" heat bath. However its possible that by 
introducing a certain frequency dependence of the coupling constants gk one can obtain 
a reasonable approximation of a heat bath on a short time scale. Such possibility should 
be subjected to further study. 

In summary, we derived the exact results for the reduced density matrix of a spin 
interacting with a type of spin bath. The precise functional dependence is determined by 
the choice of the spin bath dispersion relation and its coupling to the "central" spin. It 
turns out that the system (spin) does not reach the canonical distribution in the course 
of its evolution, and its stationary state for t = oo depends on the initial conditions. 

The author would like to thank Professors V. Privman and L.S. Schulman for their 
interest in the work and numerous fruitful! discussions. 
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